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ABSTRACT 

Deviations of the cosmic microwave background (CMB) frequency spectrum from a pure 
blackbody tell an exciting story about the thermal history of our Universe. In this paper we 
show how well future CMB measurements could decipher this tale, envisioning a PIXIE- 
like spectrometer, which could improve the distortion constraints obtained with COBE/FIRAS 
some 20 years ago by at least three orders of magnitude. This opens a large discovery space, 
offering deep insights to particle and early-universe physics, opportunities that no longer 
should be left unexplored. Specifically, we consider scenarios with annihilating and decay- 
ing relic particles, as well as signatures from the dissipation of primordial small-scale power. 
PIXIE can potentially rule out different early-universe scenarios, and moreover will allow 
unambiguous detections in many of the considered cases, as we demonstrate here. We also 
discuss slightly more futuristic experiments, with several times improved sensitivities, to high- 
light the large potential of this new window to the pre-recombination universe. 
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1 INTRODUCTION 

Cosmology is now a precise scientific discipline, with detailed the- 
oretical models that fit a wealth of very accurate measurements. 
Of the many cosmological data sets, the cosmic microwave back- 
ground (CMB) anisotropics provide the most stringent and robust 
constraints to theoretical models, allowing us to address fundamen- 
tal questions about inflation, the nature of dark matter and dark en- 
ergy, and particle physics (Smoot et al. 1992, Bennett et al. 2003, 



|Planck Collaboration et al.||2013b| >. But the CMB holds another, 
complementary piece of information: its frequency spectrum. Since 
COBE/FIRAS the average CMB spectrum is known to be extremely 
close to a perfect blackbody, with possible spectral distortions lim- 
ited to M v II V < fewx 10~ 5 JMather et al.|l994||Fixsen et al.|l996| >. 
Although thus far no CMB distortion was detected, this impressive 
measurement already places very tight constraints on the thermal 
history of our Universe, ruling out cosmologies with extended pe- 
riods of significant energy release, disturbing the equilibrium be- 
tween matter and radiation. 

More than 20 years have past since the launch of COBE and 
from the technological point of view already today it should be pos- 
sible to improve the absolute spectral sensitivity by at least three or- 
ders of magnitude (Kogut et al. 201 1). This opens a new window to 
the early Universe, on one hand allowing us to directly probe pro- 
cesses that are present within the standard cosmological paradigm, 
but also opening up a huge discovery space to unexplored non- 
standard physics. It is therefore time to ask what exactly one might 
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hope to learn from future measurements of the CMB spectrum and 
how well could the information be extracted. 

The CMB spectrum constrains energy release occurring at red- 
shifts z < z p - 2 x 10 6 . Above if,, when the Universe was only a 
few month old, the cosmological thermalization process was ex- 
tremely efficient, exponentially suppressing any distortion the ear- 
lier the energy was liberated (Zeldovich & Sunyaev|1969; Sunyaev 



& Zeldovich 1970, Illarionov & Sunyaev 1974 
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Danese & de Zotti 
1977 ; Burigana et al.|1991HHu & Silk|1993a[|Chluba|2005||Chluba 



& Sunyaev 2012; Khatri & Sunyaev 2012b). At slightly lower red- 



shift, the CMB spectrum becomes vulnerable, and any energy in- 
jection initially gives rise to a Compton y-distortion to the CMB 
blackbody, essentially by an early-universe analog of the Sunyaev- 
Zeldovich effect. At 5 x 10 4 < z < z^, the distortion rapidly evolves 
towards a chemical potential or yu-type distortion, due to efficient 
redistribution of photons over frequency, while for energy release 
occurring at z < 5 x 10 4 it keeps the shape of a -y-distortion. 

This is, however, not the end of the story. First of all, there 
is no sharp boundary at z — z M , and the more sensitive the exper- 
iment is, the deeper can one in principle probe beyond the distor- 
tion visibility function, JJ(z) - exp -(z/z,,) 5 ' 2 , defined by the ef- 
ficiency of photon production, although the challenge grows ex- 
ponentially. Secondly, the transition from p- to y-distortion oc- 
curs more gradually, with the exact shape of the distortion at the 
intermediate stages depending on the energy release history at 
10 4 < z < 3x 10 5 . Signals produced mainly during this epoch 
were discussed in Chluba & Sunyaev (2012 see Fig. 15 and 19), 
showing that the total distortion is not simply given as a super- 
position of pure p- and y-distortions. More recently, this was also 
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demonstrated by Khatri & Sunyaev (2012a) and Ch luba| ( [2013fr . 
The small (^ 10% - 30%) residuals provide additional leverage 
for distinguishing different energy release scenarios in the future, 
although no more than = 16% of the liberated energy is carried 
by this signal, making it a correction to the main mix of fi- and 
y-distortion (Chluba|2013) . 

A large number of astrophysical or cosmological processes at 
z S Z/j exist, leading to predictions of observable distortions: 

• Reionization and structure formation: the first sources of ra- 
diation during reionization (Hu et al. 1994b), supernova feedback 
(Oh et al. 2003) and structure formation shocks (Sun yaev & Zel-| 
|dovich|1972||Cen & Ostriker|1999||Miniati et al.|2000> heat the in- 
tergalactic medium at low redshifts (z < 10), leading to partial up- 
scattering of CMB photons, causing a Compton v-distortion. The 
distortion is expected to reach A/ v // y =* 10~ 7 - 10~ 6 and thus could 
be measured at a 100 cr using present-day technology (Kogut et al. 
|2011j >, teaching us about the average temperature of the intergalac- 
tic medium (e.g., |Zhang et al.|2004) , and promising a way to find 
the missing baryons in the local Universe, which otherwise are hard 
to observe ( |Cen & Ostriker|1999> . 

• Inflation: the Silk-damping of small-scale perturbations gives 
rise to both /u- and v-type distortions (Sunyaev & Zeldovich 1970. 



|Daly|1991||Barrow & Coles|1991||Hu et al.|1994a| >, which directly 
depend on the shape and amplitude of the primordial power spec- 
trum at scales 0.1 kpc < X < 1 Mpc (Chl uba et al.|20 12b; Khatri 
|et al.|2012b[ >. This allows constraining the trajectory of the infiaton 
at stages unexplored by CMB anisotropies and other ongoing or 
planned experiments ( Chluba et a l.|2012a| >. The distortion in princi- 
ple is also sensitive to the difference between adiabatic and isocur- 
vature perturbations (Barrow & Coles 1991 . Hu & Sugiyama 1994, 
Dent et al. 2012, Chluba & Grin 2013), as well as primordial non- 
Gaussianity in the squcczcd-limit. leading to a spatially varying 
spectral signal that correlates with CMB temperature anisotropies 
as large scales (Pajer & Zaldarriaga 2012; Ganc & Komatsu 2012). 

• Decaying or annihilating relics: measurements of the CMB 
spectrum also have the potential to constrain decaying and annihi- 
lating particles in the pre-recombination epoch (Hu & Silk 1993b, 
|McDonald et~aT1[200T] |Chluba||20T0} |Chluba & Sunyaev||2012f 
This is especially interesting for decaying particles with lifetimes 
t x — 2.6x 10 8 sec-2.2x 10" sec, because the shape of the distortion 
encodes when it decayed (Chluba & Sunyaev 2012). 

• Cosmological recombination radiation: the cosmological re- 
combination process of hydrogen and helium introduces distortions 
at high redshifts (z a 10 3 - 10 4 ), corresponding to ^ 260 kyr (Hi), 

130kyr (Hei), and =s 18 kyr (Heii) after the big bang ("Rubino- 



Martin et al. 2006, Chluba & Sunyaev 2006, Rubino-Martin et al. 
2008 1. The overall distortions is very small (AI V /I V ^ 10~ 9 ) but it 



has a unique frequency-dependence which opens an independent 
path to determination of cosmological parameters (like the baryon 
density and pre- stellar helium abundance) and direct measurements 
of the recombination dynamics, probing the Universe at stages well 
before the last scattering surface (Sunyaev & Chluba 2009). 

• Cooling of matter: the adiabatic cooling of ordinary matter 
continuously extracts energy from the CMB photon bath lead- 
ing to a small but indisputable distortion that directly depends on 
the baryon density and is characterized by a negative p.- and y- 
parameter at the level of - few x 10~ 9 ([Chluba 2005] |Chluba & 
|Sunyaev|20i2)|Khatri et al.|2012a) . 

All these examples demonstrate that the CMB spectrum provides 
a rich and unique source of complementary information about the 



early Universe, with the certainty for a detection of spectral dis- 
tortions at a level within reach of present day and future instru- 
mentation. The CMB spectrum could also allow placing interesting 
constraints on the power spectrum of small-scale magnetic fields 
(Jedamzik et al. 2000), primordial black holes (|Carr et aL|2010), 



cosmic strings ( Ostriker & Thompson 1987;Tashiroet al. 2012b a), 
and other new physics (Lochan et al. 2012; Bull & Kamionkowski 
|2013[ >, to mention a few more exotic examples. Deciphering all 
these signals will be a big challenge for the future, but it holds 
the potential for new discoveries, providing additional, independent 
constraints on unexplored processes that otherwise might remain a 
secret of our Universe. 

In this paper we investigate how well future measurement of 
the CMB spectrum might be able to constrain and distinguish the 
above mechanisms. In addition to the huge y-distortion introduced 
at low redshifts, we consider scenarios of decaying and annihilating 
relic particles and the dissipation of small-scale acoustics modes 
as representative examples. These mechanisms can all potentially 
lead to large distortions, compatible with existing limits, and thus 
show the potential of future CMB spectrum measurements. We fix 
the background cosmology and only vary parameters related to the 
distortions. This also means that the signal caused by the adiabatic 
cooling of ordinary matter can be predicted with very high preci- 
sion, and thus is taken out. The signatures from the recombination 
epoch are close to the detection limit of a PIXIE-type experiment, 
and we leave a more in depth discussion to some future work. 

Our projections demonstrate that PIXIE can potentially rule 
out different energy release scenarios, and will also allow unam- 
biguous detections for many of the considered examples. However, 
the constraints remain model-dependent and degeneracies between 
different scenarios exist. We also discuss slightly more futuristic 
experiments, with several-times improved sensitivities, highlight- 
ing the large discovery potential of this unexplored window to the 
early Universe 



2 MOCK SPECTRAL DISTORTION DATA 

Energy release in the early Universe causes spectral distortions 
of the CMB, with the shape of the distortion depending on how 
strongly the thermal history is affected. Assuming a specific sce- 
narios one can accurately predict the spectral distortion at differ- 
ent frequencies today (e.g., see |Chluba & Sunyaev 2012, Chluba 
|2013| >. To answer how well different scenarios can be distinguished, 
we have to produce mock spectral distortion data. We envision a 
PIXIE-like experiment, with many equidistant channels over a wide 
range of frequencies (30 GHz < v < 6THz). As a first step we 
shall assume that the measurement is only limited by uncorrelated 
instrumental noise. Foregrounds due to dust and synchrotron emis- 
sion were removed to a level below this sensitivity, making use of 
high frequency channels (v > 1 THz) and spatial templates [e.g., 
obtained with Planck (Planck Collaboration et al. 2013a)]. It seems 
that these requirements can be achieved in the future, with realistic 
error bars AI V ^ 5 x 10 ~ 26 Wnr 2 s" 1 Hz" 1 sr~' per Av ^ 15 GH z 
channel over a(^ 0.5 yr measurement periocrll Kogut et al. 201 1 1. 
The spectral distortion signal from early energy release is 
mainly important in the 30 GHz < v < 1 THz channels. Unless 
stated otherwise, we assume that all channels in this frequency 



1 Currently, PIXIE is meant to devote only as 25% of the 2 yrs observing 
time to measurements of the CMB spectrum (Kogut, priv. com.). 
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range can be used to constrain the thermal history. We furthermore 
assume that the channels are independent and described by a top- 
hat filter, VK,(v), centered at frequency v,. 

To accelerate the parameter estimation process we compute 
the distortions for different scenarios using the Green's function of 
the cosmological fhermalization problem (Chluba 2013): 



A/ V (z = 0) = I G th (v, z', 0) ry dz'. 



dz' 



(1) 



Here d(Q/p y )/ dz is the effective heating rate and p 7 oc r 4 the en- 
ergy density of the undistorted CMB. We also assume that the back- 
ground cosmology is fixed to Y p = 0.24, O m = 0.26, Q.b = 0.044, 
D A = 0.74, n k = 0, h = 0.71, and N tff = 3.046. Since we only need 
the average signal in specific band we first compute the averages of 
G,-,th(z, 0) = f Wj(v) G t h(v, z, 0) dv in channel i. This further reduces 
the total computational burden. 

To compute constraints on the different scenarios we use a 
Markov-Chain-Monte-Carlo method (MCMC), adapting routines 
that were developed as part of SZpack ( Chluba et al. 20 1 2c 20 1 3 1 
and are based on the Python packages of Foreman-Mackey et al. 
(2012). Alternatively, one can use simple Fisher forecasts, however, 
the parameter space becomes very large and degenerate close to 
the detection limit, so that we chose to follow an MCMC approach 
throughout. We place very conservative priors on the parameters 
and usually ran chains with ^ 10 s samples. Thanks to the adopted 
Green's function approach this is possible on a standard laptop in 
only a few seconds to minutes. The developed tools are now part of 
CosMoTHEtuvfl lChluba & Sunyaev|2012| >. 



2.1 Late time y-distortion 

The dominant distortion signals is caused by the heating of matter 
during the reionization epoch, introducing a y-distortion, 



2hv 3 2hv 3 

AI v /y = —p-Y^ix) = 



[xcoth(x/2)-4], (2) 



c 2 (e* - l) 2 

with Compton parameter, y as few x 10~ 7 Hu et al.|l994b i. We 
shall use y lc = 4 x 10~ 7 as fiducial value, assuming that y =s 2 x 
10~ 7 is caused by Compton scattering of CMB photons by thermal 
electrons, while another y ^ 2 x 10~ 7 is related to other heating 
mechanisms mentioned above. This is rather conservative, since the 
low redshift signal in principle can exceed the level of y ^ 10~ 6 qOh] 
|et al.|2 003 ). This will make a detection only easier and in terms of 
the error budget (as foreground for the smaller primordial signal) 
our numbers should produce pretty reliable estimates. 



2.2 Parametrization of the energy release mechanisms 

The following parametrization for the energy release mechanisms 
shall be considered (see Chluba & Sunyaev 2012 for more details): 

o decaying particle: ^^ = ^0 } e" r *< oc r 4 e -<~x /:)2 

o s-wave annihilation: — -p 2 - = H - / ann oc z 

o p-wave annihilation: — ^- = "^ — f mn oc const 

o dissipation acoustic modes: — -p 1 - °c z -1 (for scale-invariant 
primordial power spectrum of curvature perturbations) 



where N a — 1.9 x 10~ 7 (1 + z) 3 denotes the number density of hy- 
drogen nuclei and H (z) the Hubble parameter. The model param- 
eters for the first three mechanisms are the heating efficiency of 
decaying particles, f x (which is related to the mass of the relic par- 
ticle and the decay channels), the lifetime of the decaying particle, 
1 and the annihilation efficiency, / a 



tx 



(see 



Chluba & Sun- 



|yaev|2012| for more explaination). The annihilation efficiency has 
different definitions for the s- and p-wave cases (a factor kTo/nixc 2 
was absorbed for the p-wave scenarios). Furthermore, we have to 
multiply the total heating rate given above by the branching ra- 
tio, gh(z), for which we followlChen & Kamionkowski (2004b and 
Chluba et aL](|2010^ At z > 10 4 one has g h (z) S 1, while after re- 



combination gh(z) <k 1. The exact shape depends on the recombi- 
nation history and how efficiently the decay products transfer their 
energy to the medium. More sophisticated calculations for the heat- 
ing efficiencies based on different particle models can be found in, 
e.g., |Slatyer et al.H2009) and |Valdes et al.| ( |2010) , but for the pur- 
pose of this paper the approximation mentioned above will suffice. 
To describe the energy release caused by the dissipation of 
acoustic modes at small scales we follow Chluba & Grin ( j2013| >, 
including only adiabatic modes with 



d{QI Py ) 
dz 



'L 



2D 1 I Pi(k)d z e- 2k2/k idlak, 



(3) 



where P^k) = A, ■(k/k )" s - l+ 2"™ [nm ° ) is the small-scale power 
spectrum of curvature perturbations, k^(z) the dissipation scales, 
and D 1 a 0.81 the heating efficiency for adiabatic modes (assum- 
ing A'cft = 3.046). The distortion does depend on the type of initial 
conditions (adiabatic versus isocurvature), however, as shown by 
Chluba & Grin (2013), constraints can only be derived in a model- 
dependent way and thus discussion of adiabatic modes sweeps the 
whole possible parameter space (modulo overall efficiency factors 
and changes in the shape of the primordial power spectrum to ac- 
commodate for the differences between the perturbation modes). 

As cutoff scales we choose fc cut a 0. 12Mpc~', which repro- 
duces the heating rate caused by mode dissipation pretty well, even 
around the recombination epoch. To minimize the time spent on nu- 
merical integration, given the power spectrum parameters we tabu- 
late the heating rate prior to the computation of the distortions. 

2.3 Shift in the monopole temperature 

The CMB monopole temperature is known with extraordinary ac- 
curacy, T = 2.7260 ± 0.0013 K ((Fjxsen|2009]). However, the level 
of precision that might be achievable with a PIXIE-type experiment 
will dwarf this measurement. In the thermalization calculations we 
assumed T = 2.726 K. The error that is introduced by this assump- 
tion is at most a 0.05% relative to the predicted distortion, a margin 
one can comfortably live with. We must, however, take the possible 
shift in the temperature of the reference blackbody into account. At 
first order in A = AT/T, this is just a temperature shift term, but 
even the second order correction [a y-distortion (Chluba & Sunyaev 
2004)] has to be considered, since the error in the precise values of 
To itself corresponds toy =* (5x 10~ 4 ) 2 /2 s 10~ 7 . For the parameter 
estimation problem we thus add 

2hv 3 



A7 V = 



c 

2hv 3 

A 



n?\x) 



-„p' 



(*/[! + A]) 



G(x)A[i + A] + y sz (jc): 



+ <9(A 3 ) 



(4) 



CosmoTherm is available at 



www . Chluba . de/CosmoTherm 



where n pl (x) = [e x - 1] ', G(x) = xe v /[e v - l] 2 describes a temper- 
ature shift, and 5 / sz(x) a y-distortion term with x = hv/kT . For our 
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simulations we use A = 1.2x 1CT 4 as fiducial value (this is just made 
up) and then show how well one will be able to constrain it, assum- 
ing a Gaussian prior with width as 5 x 1CT 4 around it. This is very 
conservative, since a PIXIE-type experiment could determine the 
CMB monopole temperature with lcr-precision of AT = 3 nK (see 
below). We also assume that the effects caused by the superposition 
of blackbodies related to the motion-induced CMB dipole is taken 
out. This leads to a y-distortion quadrupole with y <* 2.563 x 10~ 7 , 
and a shift of the CMB monopole temperature by Ar as 0.699 yuK 
( |Chluba & Sunyaev|2004|[20TI)[Sunyaev & Khatri|2013> . 

To accelerate the parameter estimation, we again first compute 
the averages over the frequency filters. We also use this procedure 
for the annihilation and decaying particle scenarios, since the pa- 
rameter dependence is sufficiently simple. For the energy release 
caused by the dissipation of acoustic modes we explicitly integrate 
the Green's function, but we tabulate the energy release history 
once the power spectrum parameters are chosen. 



3 ANNIHILATING PARTICLE SCENARIOS 

As first simple scenario we consider an annihilating particle with 
p-wave annihilation cross section [(crv) cc T (see McDonald et al. 
|2001| for related discussion)]. Constraints on this case can be de- 
rived from BBN (due to the sensitivity of the light element abun- 
dances on the baryon to photon ratio, rj), implying that the total 
amount of energy release at that epoch cannot exceed Ap y /p y ^5% 



v - distortion 



Ll-v trnnsmon 



it - distortion 



1 Steigman 2007 1. This places a bound / ann < 4 x 10 eV sec on 



the annihilation efficiency. COBE/FIRAS constraints are a factor 
of a 3 more stringent, implying / ann < 1.5 x 10 -24 eVsec~' from 
\fi\ < 9 x 10~ 5 . The tightest constraint derives from measurement 
of the CMB temperature and polarization anisotropies, correspond- 
ing to / ann < 10~ 26 eVsec~', as we argue below. Still, this suggests 
that in principle large energy release can be accommodated for this 
scenario, without violating existing constraints. Due to the redshift- 
dependence of the heating rate, most energy is liberated during the 
;u-era (see Fig.fTJ, and hence the distortion should be easily distin- 
guishable from the large y-distortion created low redshifts. 

In Fig. [2] we show the projected constraints for a PIXIE-type 
experiment in a large p-wave annihilation cross section scenario, 
with Apylpy ez 6.2 x 10~ 7 going into the distortion. For y rc we as- 
sumed a flat prior over the interval y re 6 [0, 1.5 x 10~ s ], while we 
sampled / ann uniformly between and 100 times the input value. 
Although this prior was rather wide, the MCMC computation con- 
verged very rapidly, using about 10 5 samples. A PIXIE-type ex- 
periment will easily distinguish the associated distortion from the 
reionization signal, measuring / ann with as 1.6% (lcr-error) preci- 
sion. Since the signal is directly proportional to / ann we find 



A/„ 



/an 



as 1.6% 



./an 



10- 26 eV sec- 



A/ v 



A/ PIXIE 



(5) 



for the error, where A/f IXIE =* 5 x 10~ 26 W irr 2 s" 1 Hz 1 sr"" 1 denotes 
PIXIE's sensitivity (we confirmed this statement numerically). The 
lo--detection limit of PIXIE is roughly at /: mn , p s 3x 10" 28 eV sec -1 . 
Increasing the sensitivity 2 or 4 times might be within reach, e.g., 
by extending the total time spent on spectral distortion measure- 
ments or by slightly improving the detectors. As our results show, 
this would further tighten possible limits on this scenario, allowing 
to constraint Majorana particles annihilating into lighter fermions 
i ]Goldberg|1983) . 

Figureplalso shows that the monopole temperature and reion- 
ization y-parameter could be measured with impressive accuracy, 




v [GHz] 



Figure 1. Different s- and p-wave annihilation scenario discussed in Sect.R] 
The upper panel shows the energy release rate for all cases, while the lower 
panel only illustrates the spectral signal for the small distortion scenario. 
For comparison, we show a y-distortion of y = 2 x 10~ 9 , which for PIXIE 
sensitivity, A/ v = 5xl0~ 26 WirT 2 s Hz 1 sr" 1 , should be detectable at the 
lcr-level. An unambiguous detection of the signal from the small distortion 
scenarios will be challenging even at ^ 4 times the sensitivity of PIXIE, but 
the large distortion scenarios can be tightly constrained. The amplitude of 
the distortion signal is directly proportional to the annihilation efficiency, 
while the shape just depends on the temperature-dependence of the annihi- 
lation cross section (s-wave versus p-wave) 



corresponding to AT as 3 nK and Ay re /y rc < 1%. Both A and y re are 
anti-correlated with / ann : although the annihilation distortion sig- 
nal does not include any pure temperature shift contribution it is 
not fully orthogonal to the signal related to A [see. Eq. |4}]. Simi- 
larly, every annihilation is associated with some late energy release 
(z £j 10 4 ), during the y-era, and thus boosted annihilation efficiency 
leaves less room for contribution to y from after recombination and 
during reionization, explaining the behavior. 

Assuming a relic particle with / ann . p - 10~ 28 eV sec~', we find 
that for P/X/Zs's-sensitivity the signal is below the detection limit, 
and even at 4 times increased sensitivity, only a marginal detection 
of the distortion caused by the annihilation energy release is possi- 
ble. The measurements of A and y rc are not severely compromised 
by adding this possibility to the parameter estimation problem, be- 
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-6-1-2 2 4 6 



A* [1(T 9 ] 
(A* = A - A f ) 





Fiducial values: 
A f = 1.2 X 1(T 4 
y rc = 4xl0- 7 
/a„„ = Kr 26 eV sec" 1 



3.97 3 98 3 99 4.00 4.01 4.02 4.03 






0.94 0.96 0. 



1.02 1.04 1.06 



Am, [10" 26 eVsec _1 ] 



Padmanabhan & Finkbeiner 2005, Zhang et al. 20061, where the 
best observational limit is obtained from WMAP (Gall i et al.|2009[ 
|Hutsi et al.|2009||Slatyer et al.|2009t|Hutsi et al.|20"TT) , translating 



into / anns < 2 x 10~ 23 eV sec' 1 ( jChluba et al.|2010) . This case is as- 
sociated with an energy release of Ap r /p 7 =; 8.3xl0~ 9 , available for 
spectral distortions. In contrast to the p-wave annihilation scenario, 
energy is liberated more evenly per logarithmic redshift interval, 
so that the associated spectral distortion lies between a p and y- 
distortion (see Fig.[Tp. Annihilations with f dans a 2x 10~ 23 eV sec~' 
remain undetectable, even with 4 times the sensitivity of PIXIE, in 
agreement with conclusion from previous analyses (Chlub a~et al.| 
2010; Chluba & Sunyaev 2012). A 3<r-detection becomes possible 
with 10 times the sensitivity of PIXIE. 



On the other hand, assuming f m 



10~ 22 eV sec~ 



6cr- 



detection would be possible at 4 times PIXIE sensitivity, although 
this scenario is already in tension with CMB anisotropy constraints. 
The error for the s-wave annihilation scenario roughly scales as 



A/ a „n. S 

/ann,s 



as 17% 



./an 



10- 22 eV sec-' 



A/,. 



A/ PIXIE 



(6) 



Figure 2. Large p-wave annihilation scenario. The solid black lines show 
the constraint for PIXIE sensitivity, while the other curves are for 2 and 4 
times higher sensitivity. The contours show 50% and 95% confidence levels. 
The shaded regions are for PIXIE sensitivity only. 




A* [10 "} 

(A* = A - A f ) 





ft. [in- 



Fiducial values: 
A f = 1.2 X 10~ 4 
y re =4x 10~ 7 
/ ann , s = 10- 22 eVsec- 1 
/amu>= lO-^eVsec- 1 






f [lO^-^eV »T 






/^pIIO^'oVscc- 1 ] 



Figure 3. Large distortion s- and p-wave annihilation scenario. Contours 
and lines are as before. Degeneracies between the parameters prevent a dis- 
tinction of the signatures of both particles, even for high sensitivity. 



cause the additional signal is very small. To obtained an unambigu- 
ous 5cr-detection of the p-wave annihilation signal in this scenario, 
the sensitivity needs to be increased ^10 times over PIXIE. 

Next we consider energy release due to s-wave annihilation, 
for instance associated with a dark matter particle. The annihilation 
efficiency is already tightly constrained by the effect on the CMB 
anisotropies (Peebles et~aL||2000| [Chen & Kamionkowski||2004| 



The current limit on / ann s derived from CMB anisotropies may be 
improved by another factor of ^ 6 (e.g., see Hutsi et al"||2009| 
|2011| for projections) with the next release of Planck (which will 
include all the temperature an polarization data), ACTpol and SPT- 
pol JNiemack et al.|2010||M"cMahon et al.|2009) . At this level of 
sensitivity it will be hard to directly compete using spectral distor- 
tion measurements, however, the spectral distortion constraints are 
independent and probe different epochs of the evolution, providing 
another important handle on possible systematics, e.g., related to 
possible uncertainties in the cosmological recombination process 
(Farhang et al.||20L2| |2013). Additional freedom could be added 
due to Sommerfeld enhancement of the annihilation cross section 
(e.g., see Hannestad & Tram 2011), but a more detailed investiga- 
tion of this aspect is beyond the scope of this work. 

FigurefTlalso indicates that in the p-wave annihilation scenario 



with f m 



10 eV sec a similar amount of energy is deposited 



during hydrogen recombination (z - 10 3 ) as in the well constrained 
s-wave annihilation scenario with / anns = 2 x 10~ 23 eVsec -1 . We 
thus did not consider cases with larger p-wave annihilation cross 
section, because these would already be in tension with the CMB 
anisotropy data. Improving the limit on p-wave annihilation sce- 
narios with CMB anisotropy measurements will, however, be very 
hard and the distortion signal has a larger leverage, offering a way 
to detect the signatures from particles with p-wave annihilation ef- 
ficiency /a nn .p > few x 10~ 28 eV sec -1 at PIXIE' s sensitivity. 

Finally, in Fig. 15] we show the large distortion scenario of 
Fig. [T] with simultaneous energy release due to particles with s- 
and p-wave annihilation. The parameters becomes rather degener- 
ate, and a separate detection of the s-wave annihilation effect re- 
mains challenging even at 4 times the sensitivity of PIXIE. Con- 
sidering a scenario with more comparable distortion contributions 
from s- and p-wave annihilations (/ arm ,s - 2 x 10~ 23 eVsec~' and 



Jan 



10 eVsec ), we find that an improvement of the sensi- 



tivity by a factor of ^ 40 is needed to start distinguishing the signals 
from both particles in this case, rendering an analysis along these 
lines more futuristic. This is because for this scenarios the signal 
is close to the detection limit of PIXIE, and the differences with 
respect to a pure superposition of p- and y-distortion, which could 
be used to distinguish the two cases, are only a small correction, 
necessitating this large improvement of the sensitivity. 
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Figure 4. Large and small distortion decaying particle scenario. Contours 
and lines are as before. For large energy release the distortion can be easily 
constrained, however, for small energy release the parameter space becomes 
more complicated and higher sensitivity improves matters significantly. 



4 DECAYING PARTICLE SCENARIOS 

Decaying relic particle with lifetimes = 380 kyr (corresponding to 
the time of recombination) are again tightly constrained by mea- 
surement of the CMB anisotropies (Zhang et al. 2007, Giesen et al. 
2012), while particles with minute lifetimes can affect the light el- 
ement abundances and bounds derived from BBN apply (Kawasaki 
et al. 2005 , Jedamzik 2008). However, experimental constraints for 
particles with lifetimes ^ 10 6 - 10 12 sec are less stringent, still leav- 
ing rather large room for extra energy release Ap r /p r < 10 -6 - 10~ 5 
(e.g., |Hu&Silk|1993b||Kogutetal.|20li) . A PIXIE-type CMB ex- 
periment thus has a large potential to discover the signature of some 



long-lived relic particle, or at least provide complementary and in- 
dependent constraints to these scenarios. If most of the energy is 
released at z > 3 x 10 5 a pure p-distortion is created, so that this 
case is practically degenerate, e.g., with scenarios that include an 
annihilating particle with p-wave annihilation cross section. How- 
ever, for energy release around z — 5 x 10 4 the distortion can differ 
sufficiently to become distinguishable. 

In Fig. H] we show the projected constraints for a large and 
small distortion scenario, with energy release Ap y /p r =^ 6.4 x 10~ 6 
and Apy/p y ^ 1.3x 10~ 7 , respectively. Since the total energy release 
scales as Ap y /p r oc f x /zx (cf. |Chluba & Sunyaev|2012} , it is best 
to consider the variables fx/zx and zx — 4.8 x 10 y T x ' sec 1 ' 2 as 
parameters. This reduces the parameter covariance significantly. To 
accelerate the computation we furthermore tabulate the distortion 
for different particle lifetimes and interpolate on this grid to obtain 
the resulting distortion. With this method =* 10 5 samples can be 
taken in a few minutes on a standard quad-core laptop. 

One can see that for the large distortion scenario a - 1% preci- 
sion can be achieved for fxlzx and Zx assuming PIXIE sensitivity. 
The uncertainty on y rc increases to about ^ 3%, due to correlations 
with the signal induced by the decaying particle. This is simply 
because a noticeable fraction of the decay energy goes into produc- 
tion of y-distortions at late time, which induces an anti-correlation 
of fx/zx and y re , but a correlation of zx and y rc (increasing zx 
means less energy-release at low redshift close to recombination, 
and hence more of the y-distortion is attributed to y re ). Considering 
the small distortion scenario (with Ap y /p r ^ 1.3 x 10~ 7 going into 
distortions), shows that at PIXIE sensitivity the parameter space 
becomes rather large, showing extended regions of low probabil- 
ity due to degeneracies and correlations. Improving the sensitivity 
four times significantly tightens possible constraints on these sce- 
narios, allowing better than 5<x-detections of the particle signature. 
The constraint on the amount of energy that is released (oc fx/zx) 
is less prone to changes in the sensitivity than zx, being a proxy for 
the particle's lifetime. This is because sensitivity to zx is introduced 
mainly by the ability to distinguish a superposition of pure p- and 
y-distortion from full distortion (Chluba & Sunyaev 2012, Khatri 
& Sunyaev 2012a, Chluba 2013), but the residuals are a correction 
and thus harder to utilize. 

One can ask the question about how strongly the errors change 
when leaving the total energy release constant, but varying the par- 
ticle lifetime. For zx S= 10 4 one expects degeneracy with the y- 
distortion created at low redshifts, while for lifetimes shorter than 
t x — 3x 10 8 sec the signal becomes maximally orthogonal (p versus 
y-distortion). In Fig.Blwe illustrate this dependence of the spectral 
distortion on zx- Decreasing zx (i- e -> increasing the lifetime) moves 
the distortion from p- to a y-distortion. The residual of the distor- 
tion with respect to a superposition of pure p- and y-distortion is 
largest for zx — 8 x 10 4 , reaching roughly 30% of the total signal 
at 100 GHz. Both closer to zx - 10 4 and zx - 3 x 10 5 the residual 
becomes smaller, making a distinction harder. 

Assuming fx/zx = 1 eV and PIXIE sensitivity we find that 
for z x - 4 x 10 4 - 2 x 10 5 the errors on fx/zx and zx are typi- 
cally better than < 30%. At zx — 10 5 we obtain a =* 8% error on 
fx/zx and - 6% error on zx> representing one of the best cases. For 
Zx Si 4 x 10 4 the degeneracy with y re already becomes too large, 
and the error on zx inflates to > 27%. Similarly, for zx <: 2 X 10 5 
the signal is already too close to a pure p-distortion, which causes 
a large degeneracy between fx/zx and zx (with multi-modal so- 
lutions), simply because simultaneously increasing fx/zx and zx 
(to compensate for the suppression of the distortion amplitude by 
thermalization) gives rise to the same distortion. In other words, 
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Figure 5. Lifetime effect for different decaying particle scenarios. The up- 
per panel shows the energy release rate for all cases, while the central panel 
illustrates the distortion in comparison with a y-distortion of y = 2 X 10~ 7 . 
The lower panel shows the residual distortion after subtracting the best-fit 
fi- and y-superposition. 



a pure //-distortion is insensitive to when it was created and thus 
does not allow differentiating between scenarios with different par- 
ticle lifetimes at z > few x 10 5 . Still, a tight upper limit on the 
total amount of energy that is release can be placed, constrain- 
ing the possible abundance of decaying particles with lifetimes 
t x - 6 x 10 6 sec - 3 x 10 8 sec. 

These statements, however, depend strongly on the sensitiv- 
ity of the experiment and on how large the average distortion is. 
As explained above, the information about the particle lifetime is 
largely encoded in the deviations from a pure superposition of// and 
v-distortion, however, the residual is a correction and thus higher 
sensitivity or a larger distortion are needed to make use of that in- 
formation. Assuming fx/zx = 1 eV and zx = 2 x 10 4 , a PIXIE- 
type experiment is unable to constrain the lifetime of the particle. 
The degeneracy is already broken at twice the sensitivity of PIXIE, 
yielding =* 29% error on fx/zx and - 17% error on zx- This fur- 
ther improves to ^ 14% error on fxlzx and - 9% error on zx for 
four times the sensitivity of PIXIE. This energy release scenario 
corresponds to Ap r /p r =a 6.4 x 10~ 7 , so that the distortion is com- 
parable in amplitude to the v-signal from late times. Assuming that 
less energy is liberated by the decaying particle increases the er- 
rors (and hence the degeneracy), and conversely, for larger decay 
energy the errors diminish. Overall, a PIXIE-type experiment will 
provide a pretty good probe for long-lived particles with lifetimes 
t x - 5.8 x 10 8 sec - 1.4 x 10 10 sec and f x /z x Z 1 eV. 



5 DISSIPATION OF SMALL-SCALE ACOUSTIC MODES 

The prospect of accurate measurements of the CMB spectrum with 
a PIXIE-type experiment spurred renewed interests in how primor- 
dial perturbations at small-scales dissipate their energy (Chluba 
|& Sunyaev|20l2l |Khatri et al.|2012a||Pajer & Z aldarriaga |2012| 
Chluba et al.||2012b Dent et al.||2012l IGanc & Komatsu 2012 



Chluba et al. 2012a; Powell 2012; Khatri & Sunyaev 2013, Chluba 



& Grin 2013 I. It was shown, that this effect can be used to place 
tight limits on the amplitude and shape of the power spectrum at 
scales far smaller than what is probed with measurements of the 
CMB anisotropies, in principle allowing to discover the distortion 
signatures from several classes of early universe models (e.g., see 
|Chluba et al.|2012a| >. 

Taking a conservative perspective, one can assume that the 
power spectrum of curvature perturbations is fully determined by 
CMB anisotropy measurements at large scales, implying an ampli- 
tude A( ^ 2.2 x 10~ 9 , spectral index n$ ^ 0.96, and its running 



"rim - -0.02, at pivot scale k = 0.05 Mpc l Planck Collaboration 



et al. 2013b). This is a significant extrapolation from wavenumbers 
k < 1 Mpc -1 all the way to k ^ few x 10 4 Mpc~' , and it was already 
argued that for a PIXIE-type experiment the signal remains just 
short of the l«x-detection limit (jChluba & Sunyaev|20"T2| |Chluba| 
|et al.|20 12bi. Improving the sensitivity a few times will allow a de- 
tection of this signal, however, given that the errors on A ( , n s , and 
«mn from CMB data are now < 1%, to use spectral distortion as 
a competitive probe, factors of ^ 20 - 50 improvement are neces- 
sarjrl The strongest dependence of the distortion signal is due to 
Wrun (see Fig. [6] for illustration), since small changes affect the am- 
plitude of the small-scale power spectrum and hence the associated 
spectral distortion by a large amount (Khatri et al.|2012"a| [Chluba 



See|Powell|J2012} and|Khatri & Sunyaev|J2013| for some more in depth 



discussion of this challenge. 
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et al. 2012b I, providing some amplification. Still, this application 
of spectral distortion measurements remains futuristic, being com- 
parable to the challenge of measuring the cosmological hydrogen 
and helium recombination features with high precision. 

Both from the theoretical and observational point of view, 
there is, however, no reason to believe that the small-scale power 
spectrum is described by what is dictated by large-scale measure- 
ments. There is no shortage of models that create, bumps, kinks, 
steps, or oscillatory features in the primordial power spectrum (e.g., 



y - distortion 



11- y transition 



(l - distortion 



Salopek et al. 1989, Starobinskij 1992; Ivanov et al. 1994; Ran- 



dall et al. 1996; Stewart 1997b; Copeland et al. 1998; Starobinsky 
1998; Chung et al.||2000| |Hunt & Sarkarp007 Joy et al. 2008 






Barnaby et al. 2009 Barnaby 2010a, Ben-Dayan & Brustein 2010 



Achucarro et al. 2011, Cespedes et a l.|2012) , and direct observa 
tional constraints (e.g.. sec Bringmann ct al. 2012 for overview) 
leave large room for excess power at k > 1 Mpc~' . The recent re- 
sults obtained with Planck, e.g., from limits to non-Gaussianity 
(Planck Collaboration et al. 2013e), certainly further reduce the 
allowed parameter space for different models, but the existence 
of large-scale anomalies (Planck Collaboration et al. 20 13d I, and 
possible small-scale power spectrum features (Planck Collabora- 
|tion et al.|20 13c) indicate that matters might be more complex. A 
PIXIE-type experiment will therefore open up a new window to 
early-universe models, no matter if a distortion is detected or not. 

Given the range of possibilities, we shall pick a few illustrative 
cases, representing simple classes of models. Detailed constraints 
on specific models should be derived in a case-by-case basis, how- 
ever, our selection provides some intuition for what could be pos- 
sible in the future. We start with a simple step, AA f > 0, in the am- 
plitude of the curvature power spectrum at different k s > 1 Mpc~', 
assuming a spectral index nL. If n' s a; 1 and fc, - 1 Mpc~', from the 
practical point of view this case is degenerate with the spectral dis- 
tortion produced by s-wave annihilation [both have a heating rate 
&(Qlp y )l dz — z -1 over most redshifts]. The difference is, however, 
that a step in the small-scale power spectrum at k s =* 1 Mpc~' is 
not constrained by large-scale CMB anisotropy measurements, but 
could be tightly constrained with a PIXIE-type experiment (Chluba 
et al. 20 12a I. Similarly, degeneracy with the p-wave annihilation 
scenario exists if n' s as 2. For simplicity below we only consider the 
case n' s = n$. 

In Fig. [6] we show the heating rate and expected distortion for 
this scenario, assuming AA<- s 10~ 8 , k s = 30Mpc~', and n' s = 0.96 
for illustration. At z > 5 x 10 s the effective heating rate is about 
5 times larger than for the standard background power spectrum 
without running. Consequently, also the yu-type contribution to the 
resulting spectral distortion is found to be as 5 times larger, with 
additional contributions from the /i - y-transition era. We can also 
see that the effective heating rate changes gradually to the one of the 
background model around z - 4x 10 4 . Modes with fixed wavenum- 
ber k dissipate their energy in a range of reds hifts with a maximum 
at redshifts z diss * 4.5 x 10 5 [/t/10 3 Mpc" 1 ] 2 ' 3 jChlubaet al.|2012a i. 
Thus, no abrupt change of the heating rate is expected. 

Since the distortion in principle depends on how energy is 
released at 10 4 < z ^ 3 x 10 5 , one does expect to be sensitive 
to k s . From the discussion of decaying particle scenarios it is al- 
ready clear that only for rather large distortions (i.e., a step am- 
plitude AAf > few x 10~ 7 ) will a PIXIE-type experiment be able 
to constrain the position of the step. In Fig. [7] we show the pro- 
jected constraints on this scenario, assuming that n' s = 0.96 is fixed 
with AA ( = 4 x 10~ 7 and k s = 30Mpc _1 . Both the amplitude and 
position of the step are well constrained, with Ak s /k s ^5% and 
AAAf /AAf - 1%. Increasing the sensitivity (or similarly consider- 
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Figure 6. Effective heating rate (upper panel) and associated spectral dis- 
tortion (lower panel) caused by the dissipation of small-scale acoustic 
modes in different scenarios. For reference we show a y-distortion with 
y = 1 X 10 . For the standard power spectrum we used A<_- = 2.2 X 10~ 9 
and ns = 0.96 at pivot scale ko = 0.05 Mpc~ . All but one case are without 
running. The two scenarios with a step and bend of the primordial power 
spectrum lead to rather similar distortions (modulo and overall factor), and 
thus become hard to distinguish, although each model should be detectable 
with a PIXIE-\ike experiment at more than 5<r-confidence. 



ing a scenario with larger step amplitude) further tightens the con- 
straints. Similar to the discussion for the decaying particle case, 
moving k s closer to a; 3Mpc _1 (z diss a 10 4 ) or s 540Mpc~' 
fcdiss - 3 x 10 s ), the sensitivity on the position of the step is ex- 
pected to degrades. A PIXIE-type experiment will be most sensi- 
tive to a step at k s as 20Mpc~' - 50Mpc~' giving Afcjfc, < 6% 
and AAAJAA ( < 2% for KA ( > 4 x 10~ 7 . For k s > 150 Mpc" 1 and 
k s < 2 Mpc~' the error in £ s increases above a; 30%, although at the 
boundaries, the amplitude of the step can still be constrained rather 
precisely (^ 6% and ^1%, respectively). This is because the dis- 
tortion is rather large and only the information about the position 
of the step is lost at these limits. Improved sensitivity again helps 
breaking degeneracies, allowing to tighten the constraints on this 
scenario, and broadening the range over which the location of the 
step can be determined, analogous to the decaying particle case. 

As an additional example, we shall consider a primordial 
power spectrum with a change of the spectral index from n s to n b at 
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Figure 7. Large distortion scenario caused by a step in the small-scale 
power spectrum. Contours and lines are as before. A PIXIE-\ike experiment 
allows constraining scale and amplitude of a step in the power spectrum at 
k s a 20Mpc _1 - 50Mpc _1 with AA/ > 4 X 10" 7 to < 6% precision. 



some scale kt,, introducing a simple bend. This behavior could be 
expected from running mass models (e.g., Stewart 1997ab~| |Covi| 
& Lyt h|1999| >, or a small-scale isocurvature mode with blue spec- 
tral index that is completely subdominant at large scales (Chluba & 
Grin 2013 1. The shape of the distortion is determined by n b , since it 
sets theyu/y-distortion mixture, and how significant deviations from 
this simple superposition are. The value of k b just parametrizes the 
overall amplitude of the distortion, with spectral distortion mea- 
surements being insensitive to scenarios with k b > fewx 10 4 Mpc~' 
(Chluba et al. 2012a). Also, if n b < "s, this scenario will be hard 
to constrain, since already the signature from the standard power 
spectrum is rather small. For n b > n s , more energy is dissipated 
and hence a detection should be possible with a PIXIE-like experi- 
ment (see Fig.[6]for illustration). 

Estimates for the amplitude of the distortion can be computed 
from the model of the small-scale power spectrum, using Eq. (3} 
and Ap y /p y = f d(Q/py)/ dz er ( *l l i>> ~ dz. Furthermore, an approxi- 
mation of how large the effective p- and y-parameters are can be 
obtained using fc-space window function given in Chluba et al. 
(2012a) and Chluba & Grin (2013), or the simple approximations 
for the Green's function provided in |Chluba| (2013[ >. From these 
considerations it follows that the larger k b , the steeper does the 
small-scale power spectrum have to become for fixed experimental 
sensitivity to allow determination of the spectral index and bend 
location. Similarly, for fixed value of n b , sensitivity to the location 
of the bend is diminished, the larger k b becomes (the distortion be- 
comes smaller since less energy is liberated). 

For parameter estimations it is better to specify the ampli- 
tude of the small-scale power spectrum at some pivot scale in- 
stead of using k b . Setting the power spectrum amplitude, A b , at 
ko, b - 45Mpc~' keeps the total energy release roughly constant 
when changing n b . To determine A b (k b ), assuming no running of the 
background power spectrum, we can use A b = A^k" s " b k n ~" s k', 



To ensure that k b > 1 Mpc 



i K b "0 *0,b • 

(we shall not consider cases with both 



step and change of the spectral index here) we have the condition 
A b <2.0x 10-" 45""-'. 

To give some examples, for fib = 1.5 and k b = 3Mpc~'(Ab - 
7.2 X 10~ 9 ) the total energy available for creation of distortions is 
Apy/py =* 2.5 x 10~ 7 . The associated signal should be easily de- 
tectable with a PIXIE-like experiment, however, due to degenera- 
cies the underlying parameters are less constrained. We find that 
for 10 times PIXIE's sensitivity the errors are An b /n b as 10% and 
AA b /A b a 15% around the most probable solution, however, addi- 
tional solutions away from the input parameters were found, show- 
ing how challenging it is to constrain this scenario. 

Matters do not improve when assuming a larger change of 
the spectral index. For instance, for a scenario with n b = 2.0 and 
k b = 3Mpc~'(A b ez 2.8 x 10~ 8 ) the most probable solution was 
n b a 1.8, A b a 6 x 10~ 8 , y re <* 3.8 x 10~ 7 , and A* a 9 x 10~ 9 for 10 
times the sensitivity of PIXIE, illustrating the degeneracies of the 
parameter space, which can only be broken at very much higher 
experimental sensitivity. The distortion in this case is already dom- 
inated by a pure ^-distortion, which contains too little information 
for constraining the 2 model parameters, making this behavior plau- 
sible. In addition this case is degenerate with the p-wave annihila- 
tion scenario discussed above, indicating that model-independent 
constraints are hard to derive. 

We also considered some cases with simultaneous step and 
change of the spectral index. Scenarios with spectral index n' s a 1 
and large step amplitude are easiest to constrain, giving rise to a 
large overall signal and sufficient mixture of p-, y- and intermedi- 
ate distortion. Similarly, models with bumps introduced by some 
scenarios with particle production (e.g., Barnab y|2010b) could be 
directly constrained using spectral distortion, however, in this case 
degeneracies with decaying particles scenarios are expected. Given 
the large plausible parameter space we stopped our discussion of 
this problem at this stage, and look forward to more detailed in- 
vestigations including realistic estimates of foregrounds and other 
instrumental effects. 



6 ADDITIONAL ASPECTS 

In this section we mention a few caveats that might affect the calcu- 
lations carried out above at some significant level. In future work, 
these issues will have to be dealt with, however, they are beyond 
the scope of this paper. 



6.1 Cosmology-dependence of the Green's function 

For the parameter estimations we assumed that the background cos- 
mology does not affect the problem. This is true as long as effects 
(and errors) > 1% are considered. However, for scenarios with 
large disparity of the associated signals (large p-wave annihilation 
with simultaneous small s-wave annihilation signal), or when dif- 
ferences of the signal with respect to a simple superposition of pure 
p.- and y-distortion are important (e.g., when attempting to use the 
residual to learn something about the lifetime of a decaying parti- 
cle), this assumption could lead to an underestimation of the errors, 
correlations, and degeneracies. In this situation, the computations 
carried out with the Green's function method have to be extended 
to include its cosmology-dependence, and constraints should be de- 
rived simultaneously using CMB anisotropy, BAO, supernova, and 
large-scale structure data. This in principle can be easily achieved, 
however, will be left to future work. 
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6.2 Sensitivity and frequency resolution 

In the previous sections, we only varied the sensitivity of the ex- 
periment, assuming that there is no fundamental limitation down 
to what precision foregrounds can be separated. Factors of a few 
improvements over PIXIE might be within reach, e.g., by extend- 
ing the total time spent on spectral distortion measurements or by 
slightly improving the detectors, however, beyond that more so- 
phisticated modeling and optimization will certainly be required. 
In this context, one can also vary the number of frequency chan- 
nels. PIXIE is based on a Fourier-transform spectrometer (FTS) 
for which the spectral resolution is set by the dimension of the 
instrument. Therefore, even improving the frequency resolution 
by a factor of 2 might be challenging at this point and alterna- 
tive approaches may have to be considered. Adding more channels 
could generally help handling foregrounds and removal of narrow 
features [e.g., CO lines that pollute some of the Planck channels 
(Planck Collaboration et al.|2013a| >], but since the primordial dis- 
tortion signal is rather broad this does not lead to any obvious lever- 
age with respect to the distortion signal itself. Sensitivity appears 
to be more crucial. 

Another aspect is the lowest frequency channel. Residuals of 
the distortion with respect to a superposition of pure yu- and y- 
distortion are noticeable even below the 30 GHz channel of PIXIE 
(see Fig.Bl, so that adding bands might further help to increase the 
sensitivity to different energy release scenarios. In an FTS approach 
this is again limited by the dimension of the instrument. Ground- 
based observations might provide and interesting alternative avenue 
at those frequencies. 



6.3 Pre-recombination emission 



It was shown (Chluba & Sunyaev 2009) that the traces of neutral 
hydrogen and helium atoms present in the pre-recombination epoch 
respond strongly if the CMB spectrum is distorted. In this situation, 
cycles of uncompensated free-bound and bound-bound transitions 
lead to a significant enhancement of the cosmological recombina- 
tion radiation, giving rise to additional spectral features that can 
be much more prominent than the lines produced during the re- 
spective recombination epochs. These features might allow to tell 
the difference between post-recombination and pre-recombination 
y-distortions, thus providing another way to determine the time- 
dependence of energy-release processes at 10 3 < z < 10 4 . 

Both at low and very high frequencies these features can be- 
come as large as the primordial spectral distortion itself (see Fig. 
(12) of Chluba & Sunyaev 20091. This therefore is an important 
correction that so far has not being taken into account consis- 
tently, but has the potential to alter the thermalization efficiency at 
z ;S 2x 10 6 . In particular, yu-distortions cause strong induced atomic 
emission and thus enhances the atomic photon production rate at 
low frequencies significantly. Collisional processes (non-radiative) 
will diminish the efficiency of this process, but overall it could af- 
fect the detailed shape of the distortion caused by energy release at 
an important level. Especially, if the primordial distortion is close 
to the COBE/EIRAS limits, this problem will have to be solved if 
in the future high-precision spectral measurements will be used to 
learn details about the thermal history of our Universe. 



7 CONCLUSIONS 

We demonstrated that a PIXIE-type experiment will open an unex- 
plored window to the early Universe. It will not only be able to rule 
out different scenarios with annihilating and decaying particles, or 
large small-scale power, but could furthermore allow constraining 
specific model-parameters on a case-by-case basis. Our analysis 
shows that improving existing limits on s-wave annihilation sce- 
narios using spectral distortions will be challenging, but for p-wave 
annihilation accurate measurements (percent precision) of the an- 
nihilation efficiency will be possible if / ann . p > 10~ 26 eV sec -1 as- 
suming PIXIE's specifications (see Sect.pl. A p-wave annihilation 
signature could be detectable down to /ann.p - fewx 10~ 28 eV sec -1 , 
a sensitivity that will be very hard to achieve using, e.g., CMB 
anisotropies. Directly distinguishing the signatures from s- and p- 
wave annihilations, however, requires =s 40 times improvement 
of the experimental sensitivity over PIXIE, rendering an analysis 
along these lines more futuristic. 

A PIXIE-type experiment will also be an exquisite probe for 
long-lived particles with lifetimes f x - 5.8 x 10 8 sec- 1.4 x 10 10 sec 
and decay efficiency, fx/zx ^ 1 eV. In this case distortions in prin- 
ciple can be used to constrain both the lifetime and decay efficiency 
with percent precision. For particles with longer/shorter lifetime, 
the distortion signal is too close to a pure y/yu-distortion to allow 
telling the particles lifetime. In both cases, however, tight limits on 
the abundance and mass of the particle can be derived from spectral 
distortion measurements. The extend to which different cases can 
be distinguished is a strong function of the sensitivity, so that push- 
ing factors of a few improves matters significantly, complementing 
constraints derived from the CMB isotropies (t x — 380 yrs) and the 
light element abundances {tx — few minutes). 

Finally, CMB spectral distortion can also be used to constrain 
different early Universe models, responsible for primordial per- 
turbations of the cosmic fluid at wave numbers 1 Mpc~' < k < 
few x 10 4 Mpc~'. Extrapolating the power spectrum from large 
scales, k < 1 Mpc~', all the way to k a few x 10 4 improvements 
by a factor of ^ 20 - 50 in the sensitivity of PIXIE are needed to 
deliver a probe that is competitive with current CMB anisotropy 
data, however, the parameter space in principle is wide open, and 
stringent limits on a model-by-model basis can be obtained already 
with a PIXIE-type experiment. These limits remain rather model- 
dependent even for much higher sensitivity, and multi-modal solu- 
tions are obtained. 

Simple models with a step in the small-scale power spec- 
trum can be constrained very well, if the step appears at k s =* 
20 Mpc" 1 - 50 Mpc" 1 giving Akjk s < 6% and AAA ( /AA ( < 2% for 
step amplitudes AA<- > 4 x 10~ 7 . Similar to the decaying particles 
scenarios we find that these statements are a strong function of the 
experimental sensitivity, thus factors of a few improvement provide 
a large leverage. Constraints on specific parameters of models with 
flaring small-scale power are hard to obtain, even at much higher 
sensitivity. Still, CMB spectral distortions translate into tight inte- 
gral constraints on the inflaton's trajectory. 

Given the huge possible parameter space, a PIXIE-type exper- 
iment provides a unique opportunity to probe early-universe mod- 
els and particle physics, no matter if a primordial distortion signal 
is detected and can be linked to a specific model or not. We look 
forward to improving our estimates, accounting for more realistic 
foregrounds and other experimental aspects. 
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